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Abstract 

Nuclear reactions induced by stable and/or radioactive neutron- 
rich nuclei provide the opportunity to pin down the equation of state 
of neutron-rich matter, especially the density (p) dependence of its 
isospin-dependent part, i.e., the nuclear symmetry energy E aym . A 
conservative constraint, 32(p/p ) - 7 < E Bym (p) < 32(p/p ) 11 , around 
the nuclear matter saturation density po has recently been obtained 
from the isospin diffusion data in intermediate energy heavy-ion colli- 
sions. We review this exciting result and discuss its consequences and 
implications on nuclear effective interactions, radii and cooling mech- 
anisms of neutron stars. 
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Figure 1: Left panel: density dependence of the symmetry energy S from the contin- 
uous choice Brueckner-Hartree-Fock with Rcid93 potential (circles) , self-consistent 
Green's function theory with Reid93 potential (full line), variational calculation 
with Argonne Avl4 potential (dashed line), Dirac-Brueckner-Hartree-Fock calcu- 
lation (triangles), relativistic mean-field model (squares), and effective field theory 
(dash-dotted line) [JJ. Right panel: same as left panel for the MDI interaction with 
four x parameters as used in the IBUU04 model |1U| . 



1 Nuclear equation of state and symmetry energy 
of neutron-rich matter 

Using a parabolic approximation which has been verified by all many-body 
theories to date, the equation of state (EOS) of isospin asymmetric nuclear 
matter can be written as 

E(p,8) = E(p,5 = 0) + E sym (p)5 2 + 0(5 4 ), (1) 

where 5 = (p n — p p )/(p p + p n ) is the isospin asymmetry and E sym (p) is 
the density-dependent nuclear symmetry energy. The latter is very impor- 
tant for understanding many interesting astrophysical problems [2] , the 
structure of radioactive nuclei [HI E] > an d the dynamics of heavy-ion reac- 
tions [31 El El IH1 • Unfortunately, E syra (p) is poorly known with theoretical 
predictions diverging widely at both low and high densities as shown in the 
left panel of Fig.^for some of the most widely used microscopic many-body 
theories (§] . The theoretical uncertainties are largely due to a lack of knowl- 
edge about the isospin dependence of nuclear effective interactions and the 
limitations of existing many-body techniques. On the other hand, heavy-ion 
reactions, especially those induced by radioactive beams, provide a unique 
opportunity to pin down the density dependence of nuclear symmetry en- 



ergy in terrestrial laboratories. Indeed, significant progress has recently been 
made both experimentally and theoretically in determining the symmetry 
energy at subnormal densities jlUl 1111 1121 1131 114j . Future high energy ra- 
dioactive beams will further allow one to determine the symmetry energy at 
supranormal densities. 

To extract information about the density dependence of the nuclear sym- 
metry energy from nuclear reactions induced by neutron-rich nuclei, reliable 
reaction models that take into account the isospin degree of freedom are 
needed. One such model that has been used extensively is the isospin- 
and momentum-dependent IBUU04 transport model [15] ■ In this model, 
the single-nucleon potential, which is one of the most important inputs to 
all transport models, has the following density- and momentum-dependent 
form 

U(p, 5, p, r, x) = A u (x)^ + Mx)^ +b(^X(1- x5 2 ) 
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where x, A u (x), Aj>(x), B, C TT i, a, and A are all parameters. The corre- 
sponding momentum-dependent interaction (MDI) gives an incompressibil- 
ity of Kq = 211 MeV for symmetric nuclear matter at saturation density. 
The parameter x is introduced to mimic predictions on E syra (p) by micro- 
scopic and/or phenomenological many-body theories ^H]- Shown in the 
right panel of Fig. ^ is the density dependence of the symmetry energy for 
x = —2, —1, and 1. The momentum dependence of the single-nucleon po- 
tential, given by the last term of Eq. © , leads to reduced effective masses 
for nucleons and consequently also the in-medium nucleon-nucleon (NN) 
cross sections |131 117] . In the IBUU04 transport model, both experimental 
free-space and in-medium NN cross sections can be used. 



1.1 Constraining the density dependence of the nuclear sym- 
metry energy by isospin diffusion 

Tsang et al. recently studied the degree of isospin diffusion in the reac- 
tion 124 Sn + 112 Sn by measuring ^H] 

2Xl24 Sn+ 112 Sn - Xl24 Sn+ 124 Sn - Xll2 Sn+ 112 Sn 

Ki = Y Y W 

Al24g n+ 124 Sn - All2 Sn+ 112 Sn 

where X is the average isospin asymmetry (S) of the 124 Sn-like residue. 
Shown in the left panel of Fig. [2] is a comparison of the strength of isospin 




Figure 2: Left panel: the degree of isospin diffusion as a function of K asy (p ) for free 
(filled squares) and in-medium (open squares) nucleon-nucleon cross sections 13 . 
Right panel: symmetry energies obtained from 21 sets of Skyrme interactions and 
the MDI interaction with x = — 1 and x = |14|. 

diffusion 1 — R; obtained with free (filled squares) and in-medium (open 
squares) NN cross sections as a function of the asymmetric part of the iso- 
baric incompressibility of nuclear matter at saturation density po |1U| . i.e., 

Kasyipo) = 9pl (d 2 E sym /dp 2 y ^ - 18p (dE sym /dp) pQ . (4) 

It is interesting to note that with in-medium NN cross sections the strength 
of isospin diffusion 1 — R; decreases monotonically with decreasing value 
of x (i.e., increasing stiffness). With free-space NN cross sections, there 
appears, however, a minimum at around x = — 1. The value of K asy (pQ) 
determined from the use of in-medium NN cross sections in the IBUU04 
model is about —500 ± 50 MeV. This value is more precise than previous 
estimates |19| and thus represents a significant advance in resolving the 
major remaining uncertainty about the incompressibility of nuclear matter 
at normal density |20l I21j . Also shown in Fig. [21 are the values for 7 used 
in fitting the symmetry energy around po using the expression E sym (p) = 
32(p/po) 7 - From results obtained with the in-medium NN cross sections, 
the 7 parameter is constrained to be between 0.7 and 1.1. 

2 Constraining nuclear effective interactions 

The above constraints on the symmetry energy limit the nuclear effective 
interactions in nuclear matter. This can be seen by comparing them with 



the symmetry energies obtained from Skyrme effective interactions ^1] . The 
right panel of Fig. [2] displays the density dependence of E sym (p) for 21 sets 
of Skyrme interaction parameters. Most of these effective interactions are 
not compatible with that extracted from the isospin diffusion data, given by 
the MDI interaction with x parameter between —1 and 0, except the SIV, 
SV, G CT , and Gr Skyrme effective interactions. 



3 Constraining the slope of the symmetry energy 




Figure 3: Left panel: symmetric (Po), asymmetric (P a sy), and total pressure in 
neutron stars at /3-equilibrium using the MDI interaction with x — and x = — 1. 
Right panel: strength of isospin diffusion in 124 Sn+ 112 Sn reactions and size of 
neutron skin in 208 Pb as functions of the slope of the symmetry energy. The inset 
is the correlation between isospin asymmetry and density of matters at the instant 
of 20 fm/c in the considered reaction |22| . 

While most properties of neutron stars depend on both the isospin- 
independent and -dependent parts of the nuclear equation of state, their 
radii are primarily determined by the isospin asymmetric pressure that is 
proportional to the slope of the symmetry energy E' sym (p) [Q. For the sim- 
plest case of a neutron-proton-electron (npe) matter in neutron stars at (3 
equilibrium, the pressure is given by 

P(p,6) = P (p)+P asy (p,5)=p 2 ^—j^ + - Pe p e 

E'(p, 5 = 0) + E' sym (p)8 2 } + - 5) P E sym (p) (5) 



P 2 



where the electron density is p e = ^{^ ~ 8)p an d the chemical potential is 



= Mn ~ fJ-p = 45£' sym (p). The equilibrium value of 5 is determined by the 
chemical equilibrium and charge neutrality conditions, i.e., 5 = 1 — 2x p with 

X p « 0.048 [£ sy m(p)/£sym(A))] 3 (p/po)(l " 2x p ) 3 . (6) 

Because of the large <5 value in neutron stars, the electron degenerate pres- 
sure is small. Moreover, the isospin symmetric contribution to the pressure 
is also very small around normal nuclear matter density as E'(po, 5 = 0) = 0. 
Shown in the left panel of Fig. Elare the isospin symmetric (Po) an d asym- 
metric (P a sy) as well as the total pressure in neutron stars at /^-equilibrium 
calculated from Eq.© using the MDI interaction with x = and x = —1. It 
is seen that up to about 2.5po for x = — 1 and about 1.5po for x = the total 
pressure is dominated by the isospin asymmetric contribution. Because neu- 
tron star radii are determined by the pressure at moderate densities where 
the proton content of matter is small, they are very sensitive to the slope 
of the symmetry energy near and just above po- I n particular, a stiffer 
symmetry energy is expected to lead to a larger neutron star radius. 

For isospin diffusion in heavy- ion reactions, it is a re-distribution of 
isospin asymmetries that is initially carried by the colliding nuclei. The 
degree and rate of this process depend on the relative pressures of neutrons 
and protons, namely the slope of E syra (p). With a stiffer E sym (p), it is harder 
for neutrons and protons to mix, leading thus to a smaller /slower isospin 
diffusion. Because of isospin asymmetric pressure, dilute neutron-rich clouds 
surrounding a more symmetric dense region are dynamically generated in 
heavy-ion reactions through isospin diffusion as illustrated in the inset of 
the right panel of Fig. EJ where the correlation between isospin asymmetry 
and density of matters at the instant of 20 fm/c in 124 Sn+ 112 Sn reactions 
is shown. A similar effect is seen in neutron-rich nuclei, where the isospin 
asymmetric pressure pushes extra neutrons further out of an isospin sym- 
metric core of near normal density, leading to an increasing of the sizes of 
neutron-skins in neutron-rich nuclei with increasing slope of the symmetry 
energy [1|31|23]. 

From above discussions, one expects that the radii of neutron stars, the 
degree of isospin diffusion in heavy-ion collisions, and the sizes of neutron- 
skins in heavy nuclei are all correlated through the same underlying nuclear 
symmetry energy. This is demonstrated in the right panel of Fig. |H1 where 
the strength of the isospin diffusion 1 — Ri, calculated with the IBUU04 
model with in- medium NN cross sections, and the thickness of neutron skin 
dR np in 208 Pb, calculated using the Skyrme Hartree-Fock with interaction 
parameters adjusted to give an EOS which is similar to the effective in- 
teraction used in the IBUU04 model P2|, are examined simultaneously as 



functions of the slope parameter L = 3po(dE sym /dp) po of the symmetry 
energy. It is seen that 1 — R; decreases while dR np increases with increas- 
ing L as expected. Taken the fiducial value dR np = 0.2 ± 0.04 fm, that is 
measured and supported strongly by several recent calculations [2], and the 
NSCL/MSU data 1 - Ri = 0.525 ± 0.05 [IT], the L parameter is constrained 
in a common range between 62.1 MeV (x = 0) and 107.4 MeV (x = —1) 22 . 

4 Constraining the radii and cooling mechanisms 
of neutron stars 




Figure 4: Left panel: mass and proton fraction of neutron star and the symmetry 
energy as functions of density. Right panel: correlations in neutron star masses and 
radii |22|. All results are for spherically-symmetric, non-rotating, non-magnetized 
neutron stars consisting of npefj, matter at zero temperature. 

The dependence of some of the basic properties of neutron stars on the 
nuclear EOS is shown in the left panel of Fig. 0J The top and middle 
panels give, respectively, the mass of neutron star and its proton fraction Xp : 
calculated from the MDI interaction with x = 0,-1, and —2, as functions of 
its central density. The symmetry energies obtained from these interactions 
are shown in the lower panel together with that from the AV18+5v+\JTK* 
interaction of Akmal et al. (APR) [21]. It is interesting to see that up to 
about 5po the symmetry energy predicted by APR agrees very well with 
that from the MDI interaction with x = and is thus within the range of 
symmetry energies constrained by both isospin diffusion and 208 Pb neutron- 
skin thickness data. Because of the lack of experimental information on the 



EOS of pure neutron matter, the APR prediction has been widely used as 
a standard in calibrating other microscopic calculations. 

For x p below 0.14 pQ, the direct URCA process for fast cooling of proto- 
neutron stars does not proceed because energy and momentum conservation 
cannot be simultaneously satisfied. For EOSs from the MDI interaction 
with x = — 1 and x = —2, the condition for direct URCA process is fulfilled 
for nearly all neutron stars above 1 M . For the EOS from x = 0, the 
minimum density for the direct URCA process is indicated by the vertical 
dotted line, and the corresponding minimum neutron star mass is indicated 
by the horizontal dotted line. Neutron stars with masses above 1.39 M are 
thus expected to have a central density above the threshold for the direct 
URCA process. 

The relations between neutron star masses and radii for above EOSs are 
given in the right panel of Fig. 0] Also given are the constraints due to 
causality as well as the mass-radius relations from estimates of the crustal 
fraction of the moment of inertia (A//I = 0.014) in the Vela pulsar [25] 
and from the redshift measurement from Ref. [2E1- Allowed equations of 
state should lie to the right of the causality line and also cross the other 
two lines. The hatched regions are inferred limits on the radius and the 
radiation radius (the value of the radius which observed by an observer at 
infinity) defined as Roc = R/y/l — 2GM/Rc 2 for a 1.4 M neutron star. It 
is seen that the symmetry energy affects strongly the radius of a neutron 
star but only slightly its maximum mass E2j . These analyses have led to 
the conclusion that only radii between 11.5 and 13.6 km (or radiation radii 
between 14.4 and 16.3 km) are consistent with the EOSs from the MDI 
interaction with x = and x = — 1 and thus with the laboratory data on 
isospin diffusion and neutron skin thickness of heavy nuclei |22| . 

The observational determination of the neutron star radius from the 
measured spectral fluxes relies on a numerical model of the neutron star 
atmosphere and uses as inputs the composition of the atmosphere, a mea- 
surement of the distance, the column density of x-ray absorbing material, 
and the surface gravitational redshift. Since many of these quantities are dif- 
ficult to measure, only a paucity of radius measurements are available. Cur- 
rent estimates obtained from recent x-ray observations give a wide range 
of results. It is therefore very useful to compare our results with recent 
Chandra/XMM-Newton observations. Assuming a mass of 1.4 M , the in- 
ferred radiation radius Roo (in km) is 13.5 ±2.1 [2B] or 13.6 ±0.3 |29| for the 
neutron star in u Cen, 12.8 ±0.4 in M13 [3U], 14.5±^ for X7 in 47 Tuc [SI], 
and 14.5118 in M28 [32] • Except the neutron star in M13 that has a slightly 
smaller radius, all others fall into our constraints of 14.4 km < i?oo < 16.3 



km within the observational error bars that are larger than the above range 
in most cases. 

5 Summary 

The EOS of neutron-rich matter is fundamentally important for both nuclear 
physics and astrophysics. Nuclear reactions induced by neutron-rich nuclei 
provide a great opportunity to constrain the EOS of neutron-rich matter. 
The recent analysis of isospin diffusion data in heavy-ion collisions within the 
IBUU04 transport model allowed us to constrain the density dependence of 
the nuclear symmetry energy in a narrow range. Together with experimental 
information on the neutron-skin thickness of 208 Pb, this has led to stringent 
constraints on not only the nuclear effective interactions but also the cooling 
mechanisms and radii of neutron stars. 
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